Band structure effects on the interaction of charged particles with solids 
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A survey is presented of current investigations of the impact of band structure effects on various 
aspects of the interaction of charged particles with real solids. The role that interband transitions 
play in the decay mechanism of bulk plasmons is addressed, and results for plasmon linewidths in 
Al and Si are discussed. Ab initio calculations of the electronic energy loss of ions moving in Al and 
Si are also presented, within linear response theory, from a realistic description of the one-electron 
band structure and a full treatment of the dynamic electronic response of valence electrons. Both 
random and position-dependent stopping powers of valence electrons are computed. 
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I. INTRODUCTION 

A quantitative description of the interaction of charged 
particles with solids is of fundamental importance in a va- 
riety of theoretical and applied areas . When a moving 
charged particle penetrates a solid material, it may lose 
energy to the medium through various elastic and inelas- 
tic collision processes that are based on electron excita- 
tion and nuclear recoil motion in the solid. While energy 
losses due to nuclear recoil may become dominant at very 
low energies of the projectile in the case of electrons 
or ions moving with non-relativistic velocities that are 
comparable to the mean speed of electrons in the solid 
the most significant energy losses are due to the genera- 
tion of electronic excitations, such as electron-hole pairs, 
collective oscillations, i.e., plasmons, and inner-shell ex- 
citations and ionizations. 

For many years, theoretical investigations of valence- 
electron excitations in metals have been carried out on 
the basis of the so-called jellium model of the solid. 
Within this model, valence electrons are described by 
a homogeneous assembly of free electrons immersed in a 
uniform background of positive charge, the only parame- 
ter being the valence-electron density hq or, equivalently, 
the so-called electron-density parameter rg defined by the 
relation I/uq — (4/3) tt (r^ ao)"^, ao being the Bohr ra- 
dius. On the other hand, a widely used approach to treat 
long-range Coulomb interactions in a many-electron sys- 
tem is the so-called random-phase approximation (RPA) 
a mean-field theory that has been proved to be suc- 
cessful in describing both electron-hole and plasmon ex- 
citations. For instance, the free-electron gas (FEG) pre- 
dicts, within RPA, a quadratic plasmon dispersion uj{q) 
with momentum transfer which is in qualitative agree- 
ment with many experimental results Q. 

Among the limitations of a free-electron-gas descrip- 
tion of the solid there is the absence, within RPA, of 
any damping mechanism other than the Landau damp- 
ing where plasmons can decay, for momentum transfers 
above the critical value qc, through the creation of an 



electron-hole pair. As a result, the FEG predicts, within 
RPA, an infinitely sharp plasmon line for momentum 
transfers smaller than Qc, while measured plasmon lines 
of real solids show a finite linewidth Q . 

Plasmon-decay mechanisms are due, within a free- 
electron-gas description of the solid, to the existence 
of frequency-dependent many-body interactions beyond 
RPA. These high-order interactions yield non-zero prob- 
abilities for the plasmon to decay into two electron-hole 
pairs or into a plasmon of lower energy plus an electron- 
hole pair, thereby satisfying momentum and energy con- 
servation for arbitrarily small values of q. However, this 
plasmon damping, which has been shown to be of order 

for small values of the wave vector, cannot explain the 
observed finite width for q — > 0. Also, calculations by 
DuBois and Kivelson and by Hasegawa and Watabe 
Q show that this contribution to the plasmon linewidth 
is too small to account alone for the linewidth dispersion 
in simple metals. 

Additional plasmon-decay mechanisms, which yield a 
finite plasmon linewidth at zero and non-zero momen- 
tum transfers, are due to either phonon/impurity as- 
sisted electron-hole excitations or to band-structure ef- 
fects. The former can be accounted by simply introduc- 
ing a phenomenological relaxation time in the particle- 
conserving RPA- like Mermin dielectric function Q . The 
later requires the calculation of the wave-vector and fre- 
quency dependent dielectric matrix in the reduced trans- 
lational symmetry of the real solid. Early calculations 
of band-structure effects on the plasmon lifetime in sim- 
ple metals were performed by Sturm and Oliveira [p|-p^ , 
with the use of a nearly free-electron pseudopotential 
theory. Ab initio calculations of the plasmon lifetime 
have been carried out only very recently for the case of 
potassium , thereby explaining the experimentally ob- 
served anomalous dispersion of the plasmon linewidth in 
this simple metal | p2| . 

Calculations of the electronic stopping power of solids, 
i.e., the energy that moving charged particles loose per 
unit path-length due to electronic excitations in the solid. 
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have also been carried out, for many years, on the basis 
of a jellium model of the target . The calculation of 
the stopping power of real solids from the knowledge of 
the band structure and the corresponding Bloch eigen- 
functions is a laborious problem, and early theoretical 
investigations were based on either semiempirical treat- 
ments of the electronic excitations in the solid 14 -T^] 
or approximate schemes based on the assumption that 
electrons are individually bound by harmonic forces [ p^ . 
Among the most recent attempts to introduce the full 
electronic band structure in the evaluation of the elec- 
tronic stopping power for low projectile velocities there 
is, for alkaline metals, a one-band calculation |20|, as well 
as a calculation based on a linear combination of atomic 
orbitals (LCAO) j2^. The low- velocity limit was also in- 
vestigated, in the case of silicon, on the basis of a static 
treatment of the density-response of the solid |^^, and, 
more recently, ab initio band structure calculations that 
are based on a full evaluation of the dynamical density- 
response of the solid have been carried out for aluminum 
IIH and silicon ||. 

In this paper we summarize current investigations of 
the impact of band structure effects on various aspects 
of the interaction of charged particles with real solids. In 
Sec. II, we describe our full treatment of the dynami- 
cal density-response of valence electrons which is based, 
in the framework of time-dependent density-functional 
theory (TDDFT) ^j^, on a reahstic description of 
the one-electron band structure and first-principles pseu- 
dopotential theory. In Sec. Ill, band-structure effects on 
the energy-loss function are investigated. For wave vec- 
tors below the plasmon cutoff Qc we focus on both the 
dispersion of plasmon energy and plasmon linewidths, 
and for larger wave vectors we investigate the so-called 
dynamic structure factor. In Sec. IV, the electronic stop- 
ping power of valence electrons in aluminum and silicon 
is considered, which we evaluate within either random 
or channeling conditions. In Sec. V, our conclusions are 
presented. In the present work, we restrict our atten- 
tion to the lowest order in the projectile charge, thereby 
high-order corrections [p8| being ignored. Hence, existing 
differences between the stopping power for protons and 
antiprotons are not accounted here, and we focus on the 
impact of band structure effects on the stopping power 
of valence electrons in real solids. 



II. LINEAR RESPONSE 

Take a system of interacting electrons exposed to 
an external potential V^^*{r,uj). According to time- 
dependent perturbation theory and keeping only terms 
of first order in the external perturbation, the electron 
density induced in the electronic system is found to be 
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J dr'x(r,r';^)F^-*(r',c.) 
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where x(r,r';u;) represents the so-called linear density- 
response function 



X(i-,r',w) p,*o(r)p„o(i-') 



1 



Eo~En + h{uj + IT]) 
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Eo + En + h{uj + irj) 



(2) 



77 being a positive infinitesimal. /9„o(r) represent matrix 
elements, taken between the unperturbed many-particle 
ground state |^o) of energy Eq and the unperturbed 
many-particle excited state |^„) of energy En, of the 
particle-density operator 
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(3) 



where describe electron coordinates. The inverse di- 
electric function is connected with the density-response 
function by the following relation: 

e-i(r,r';cj) = (5(r-r') + J dr"u(r - r")x(r", r'; c^), 

(4) 

where i;(r — r') is the bare Coulomb interaction. 

In the framework of TDDFT, the exact density- 
response function x(r, r';w) obeys the integral equation 



X(r,r';w) = x°(r,r';t^) + J dri J drs x°(r,ri;w) 

X [w(ri -r2) + X""(ri,r2;c^)]x(r2,r';w), (5) 

where x''(r, r'; cj) represents the density-response func- 
tion for non-interacting Kohn-Sham electrons p9[| , and 
K^'^{r,r';uj) represents the reduction in the electron- 
electron interaction due to the existence of dynamical 
short-range exchange and correlation (XC) effects. 

In the case of a homogeneous electron gas, one intro- 
duces Fourier transforms and writes 

where Vq represents the Fourier transform of the bare 
Coulomb interaction, y2, , , is the well-known function of 
Lindhard ||30[], and 

~ ^VqGq^uj, (7) 

Gq.ui being the so-called local-field factor . 

For a periodic crystal, we introduce the following 
Fourier expansion of the density-response function: 

BZ 

X(r,r';c.) = ^J2I1 e'('»+«)-e-'('^+«')- Xg,G' (q, c). 
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where V represents the normahzation volume, the first 
sum runs over q vectors within the first Brillouin zone 
(BZ), and G and G' are reciprocal lattice vectors. In- 
troduction of Eq. (^) into Eqs. (^) and (||) yields ma- 
trix equations for the Fourier coefhcients xg,G'('1i ^) ^^'^ 
Eq^q, (q, w). In the case of non-interacting Kohn-Sham 
electrons, one finds: 

/ \ 1 /k.ri ^ /k+q,ri' 

Xg,G' (q> -) = 77 L L e. .„ - E.^. + h(u 



k n.n' ' ^ 



X((/)k.n e 



^i(q+G)-r| 



'k+q,n'/ Wk+q^n' 



,i(q+G')-r 



\Kn)- (9) 



The second sum runs over the band structure for each 
wave vector k in the first Brillouin zone, /k,,i repre- 
sents the Fermi-Dirac distribution function, and (/)k,Ti(r) 
and -Ek,n are Bloch eigenfunctions and eigenvalues of 
the Kohn-Sham Hamiltonian of density- functional theory 
(DFT) 1^^, which we evaluate within the local-density 
approximation (LDA) with use of the parametrization of 
Perdew and Zunger 1 34 1 . All results presented in this ar- 
ticle have been found to be well converged with the use 
in Eq. (|) of 60 (Al) and 200 (Si) bands. The sum over 
the Brillouin zone has been performed with the use of 
lOx lOx 10 and 8x8x8 Monkhorst-Pack meshes Q 
for Al and Si, respectively. 

For the evaluation of the one-electron 0k,n(r) eigen- 
functions, we first expand them in a plane- wave basis. 



0k,«(r) 



1 



G 



Wk,n(G)e- 



i(k+G)T 
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with a kinetic-energy cutoff that varies from 12 Ry in 
the case of Al (~ 100 G- vectors) to 16 Ry in the case 
of Si 300 G-vectors). Then, we evaluate the coeffi- 
cients Mk.n by solving the Kohn-Sham equation of DFT 
with a full description of the electron-ion interaction that 
is based on the use of an ab initio non-local, norm- 
conserving ionic pseudopotential ||36|| . Finally, we eval- 
uate from Eq. (^) the Fourier coefficients of the non- 
interacting density-response function, and solve a matrix 
equation for the Fourier coefficients of the interacting 
density-response function, which we obtain within either 
RPA or TDLDA. 

In the RPA, the XC kernel K^'^{r,r';LL;) entering Eq. 
(^ is set equal to zero. Within TDLDA, which represents 
an adiabatic extension to finite frequencies of the LDA, 
one writes 



dn 



(11) 



=no(r) 



where Vxc{n) represents the LDA XC potential entering 
the Kohn-Sham equation of DFT, and no(r) is the actual 
density of the electron system. In the case of a homoge- 
neous electron gas of density tiq, the Fourier transform 
of the XC kernel of Eq. (|l^) can be written in the form 
of Eq. (|) with 



where 



A 



1 Anal 
4 e'^qj. 



dVc{n) 



dn 



(12) 



(13) 



n— no 



Vc{n) being the correlation contribution to Vxc, and qp, 
the Fermi momentum. 



III. ENERGY-LOSS FUNCTION 
A. Theory 

The Hamiltonian describing the interaction between a 
particle of charge Zie at point r and the many-electron 
system is given by 



Hi = -Zxe"- \ dt 



2 / P(r') 



r - r' 



(14) 



where p(r) represents the particle-density operator of Eq. 

Within first-order perturbation theory, the probability 
per unit time for the probe particle to scatter from an 
initial state |z > of energy e.i to a final state |/ > of 
energy ey, by carrying the Fermi gas from the many- 
particle ground state | ^'o > of energy Eq to some excited 
many-particle state > of energy is given by the 
following expression: 

- y ^ |< /|i^/|*o i >|' '5(£» - £/ + So - S„) 

n 

(15) 

If one chooses the probe particle to be described by 
plane-wave states 



Mr) = -^e'""' 



and 



Mr) - ^e--, 



(16) 



(17) 



one finds the probability per unit time for the probe par- 
ticle to transfer momentum hq = Ti{q.i — q/) to the Fermi 
gas to be given by the following expression: 



P^ = j^ZfvlS{ci,^), 
where hui represents the energy transfer. 



(18) 



(19) 
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and ^(q, uj) is the so-called dynamic structure factor of 
the many-electron system: 



(20) 



Here, pq represents the Fourier transform of the particle 
density, (Pq)no are matrix elements taken between the 
many-particle eigenstates |^'o > and >, and hujno = 

En — Eq. 

The double differential cross section d^a/dfldu for in- 
elastic scattering of charged particles is simply the num- 
ber of particles being scattered per unit time, unit solid 
angle, and unit frequency into the solid angle O with 
energy transfer hiv, divided by the initial particle flux. 
From Eq. dl^), one easily finds 



dV 



dfldio An'^a, 



Zl If 



2^4 



(21) 



While high-energy transmission electron beams have 
been used to probe ^(q, uj) for q < qp, thereby providing 
experimental evidence of collective excitations Q , higher 
values of q have been studied with use of inelastic X-ray 
scattering pTf . Within the lowest-order Born approxima- 
tion, the inelastic X-ray scattering cross section is given 
by 



— — = (eo • ei)^7rr^ — 5(q, lj) 
ds Iduj ujo 



(22) 



where (eo,a;o) and (eo,wo) are the polarization and fre- 
quency of the incident and scattered photon, respectively, 
and ro is the classical electron radius. 

The dynamic structure factor, which determines the 
fluctuation of the particle density, is connected to the 
imaginary part of the density-response function of Eq. 
(^ through the fluctuation-dissipation theorem. At zero 
temperature, one writes |38|| : 

S'(q,u;) = -J j dr j dr'e-"i('-'-''lmx(r, r'; cj). (23) 

Hence, in the case of a uniform electron gas one easily 
finds the absorption probability Pq to be given by the 
following expression: 



Pr, 



hv 



Zfv^Im [-e-i(q,c^)] , 



(24) 



where Im [— e"^(q, a;)] is the so-called energy-loss func- 
tion of a homogeneous electron gas. For periodic crystals, 
introduction of Eq. (0) into Eq. (p3|) yields: 



P, 



q+G 



nv 



Ztv. 



'q+G 



Im 



^G,G 



(25) 



where q represents a wave vector in the first BZ, ?i(q-|-G) 
is the momentum transfer, and eQ^Q,(q, w) is the inverse 
dielectric matrix: 



■^G.G' 



(q, u) = Sg.g' + Vq+G XG,G' (q, w) 



(26) 



If recoil of the probe particle can be neglected, as in the 
case of heavy projectiles or fast electrons with small val- 
ues of the momentum transfer, the energy transfer is sim- 
ply huj — ?iq • V. 

We note that as long as the probe particle can be 
described by plane waves, the absorption probability 
in periodic crystals is proportional to the imaginary 
part of one diagonal element of the energy-loss matrix 

Im — eQ^Q,(q, w) . Hence, crystalline local-field effects 

[^,40|, i.e., coupHngs between q -f G and q+G' wave 
vectors with G 7^ G' only enter through the dependence 
of the diagonal elements of the inverse dielectric matrix 
on the off-diagonal elements of the direct matrix. 

When the probe particle moves with constant velocity 
V on a definite trajectory at a given impact vector b, 
initial and final states can be described in terms of plane 
waves in the direction of motion and a S function in the 
transverse direction As a result, the probability per 
unit time for the probe particle to transfer momentum ?iq 
to the Fermi gas is now given by the following expression: 

Pq = Zl ^ «q«q, ^"^^^^^^ 5(q, q' ; ^) 5,^ , 

q' 



where a; = q • v, 5„ 



(27) 

represents the Kroenecker 8 sym- 
bol, qz and are components of q and q' in the direction 
of motion, and 

S'(q,q';tj) = ^ dx ^ dr'e"' ('i''+'i'"'')lm x(r, r'; tj). 

(28) 

In the case of a homogeneous electron gas the proba- 
bility Pq of Eq. (27) is easily found to be independent of 
the impact vector and to coincide with that of Eq. (p4). 
However, for periodic crystals introduction of Eq. ^^) 
into Eq. ( p8| ) yields the position-dependent absorption 
probability: 



P, 



q+G 



2 

W 



zlY.' 



K 



^G,G+K 



(29) 



where the prime in the summation indicates that it is 
restricted to those reciprocal-lattice vectors that are per- 
pendicular to the velocity of the projectile, i.e., K • v = 
0. The most important contribution to the position- 
dependent probability of Eq. ( p^ ) is provided by the 
term K = 0, the magnitude of higher-order terms de- 
pending on the direction of the velocity. For those di- 
rections for which the condition K • v = is never satis- 
fied we have the absorption probability of Eq. (p5|), and 
for a few highly symmetric or channeling directions non- 
negligible corrections to the random result are expected. 
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We also note that the average over impact parameters of 
the position-dependent probability of Eq. ( P9[) along any 
given channel coincides with the absorption probability 
ofEq. dH). 



B. Results and discussion 

First of all, we examine the dispersion of the plasmon 
energy and linewidth by calculating the electron energy- 



loss function Im 



(q, Lo) for wave vectors that are 
For these wave vectors the 



below the plasmon cutoff q, 
FEG predicts, within either RPA or TDLDA, sharp plas- 
mon lines, which would be broadened through the intro- 
duction of a finite value of 77 in Eq. (^). In order to avoid 
numerical broadening, we have computed the energy-loss 
function of the real solid for imaginary frequencies and 
have obtained it on the real axis with 77 ^ 0+ by using 
Fade approximants [^l|j4^ ]. 

Our computed RPA plasmon peaks in Al are displayed 
in Fig. 1 , for increasing magnitude of the wave vector in 
the (100) direction. Vertical lines are 5 functions corre- 
sponding to plasmon excitation in a FEG. As a result 
of the actual band structure of the solid, the energy- 
loss function exhibits a sharp but finite plasmon peak 
at low wave vectors with a full width at half maximum 
(FWHM) that increases rapidly, for increasing q, while 
the peak height decreases according to the /-sum rule. 
For periodic crystals, one writes: 



d uj ljItci 



■ nc-G', 



(30) 



where tiq represents the Fourier components of the den- 
sity, and nG=o equals the average electron density of the 
crystal. 

RPA and TDLDA plasmon peak positions in Al, as 
calculated along the (100) direction, are plotted in Fig. 
2 by solid and dotted lines, respectively, as a function 
of the magnitude of the wave vector. The calculated 
plasmon dispersions agree with those of previous calcula- 
tions 1^,0 . For comparison, RPA and TDLDA plasmon 
dispersions corresponding to a FEG with — 2.07 are 
represented by dashed and dashed-dotted lines, respec- 
tively, showing that band-structure corrections result in a 
nearly-q-independent downward shift of plasmon energy. 
If short-range XC effects, which reduce the plasma fre- 
quency for large values of the wave vector, are included 
within the TDLDA, then one finds a result for the plas- 
mon dispersion that is in excellent agreement with ex- 
periment 

The plasmon lifetime is defined as the inverse of the 
FWHM of the energy- loss peak, A£'i/2. RPA and 
TDLDA calculations of A£'i/2 along the (100) direction 
in Al are exhibited in Fig. 3, as a function of the mag- 
nitude of the wave vector. In the limit as q — > we find 
A£^i/2(0) ~ 0.2 eV, below the early calculation carried 



out by Sturm with use of the nearly free electron ap- 
proximation [||. As for plasmon decay from the creation 
of two electron- hole pairs, which is not included in our 
band-structure calculation, one finds to lowest-order in a 
small-q expansion ]lO| ]: 



b {q/qpY ^^p 



(31) 



where b ~ 0.03 for ~ 2. Adding band-structure effects, 
as derived from our calculated RPA plasmon peaks, and 
many-body corrections, as obtained from Eq. (^l|), we 
have found the results plotted in Fig. 3 by open tri- 
angles. Sohd circles correspond to the result of not in- 
cluding plasmon decay from many-body effects, thereby 
showing that these effects give only a minor contribu- 
tion to the plasmon linewidth. As phonons are expected 
to give a relative contribution to AEi^2{q) that changes 
little with the magnitude of the wave vector , the re- 
duced plasmon linewidth AEi/2{q) / AEi/2{0) is plotted 
in Fig. 4 by open circles (RPA) and triangles (TDLDA) , 
showing a reasonable agreement with the experimental 
results of Gibbons et al [Q. 

In Fig. 5 we show full-band-structure calculations of 
RPA and TDLDA energy-loss functions in Si, as obtained 
along the (111) direction. As in Fig. 1, vertical lines rep- 
resent S functions corresponding to plasmon excitations 
in a FEG, now with ~ 2.01. Besides the finite width 
of the energy-loss function of real Si, which is found to 
be about 5 times wider than in Al, our calculated plas- 
mon energy shows a characteristic q dependence which is 
considerably different from that of a FEG. 

The energy positions of RPA and TDLDA plasmon 
peaks of Si are plotted in Fig. 6 along the (111) direc- 
tion, as a function of the magnitude of the wave vector, 
together with the FEG prediction and the experimen- 
tal result 1^^. In the long- wavelength limit {q 0), 
the plasmon energy is expected to be only slightly larger 
than in the case of a FEG, in agreement with the Penn 
model dielectric function |48| for semiconductors which 



predicts a plasmon frequency ui 



-pj ,.»-g being 
the energy of the effective band gap and LUp representing 
the plasma frequency of the corresponding FEG. We also 
note that while in the case of Al band-structure effects 
result in a nearly-q-independent downward shift of the 
FEG curve, interband transitions in Si flatten the plas- 
mon curve with a considerable reduction of the plasma 
frequency at large values of the wave vector. Though our 
calculated plasmon dispersion in Si is at low wave vectors 
in good agreement with the experiment, at large |q| our 
predictions start to deviate from the experimental value, 
a result also found in previous calculations |4S|j50[ |. 

In Fig. 7 we show our calculation of the plasmon 
linewidth dispersion of Si, as obtained along the (111) 
direction with no inclusion of crystalline local-field ef- 
fects. In the limit as g ^ we find Ai^i/jfO) ~ 1 eV, 
well below the prediction of Louie et al pl| who per- 
formed a calculation of the energy-loss function in Si at 
q = with full inclusion of local-field effects. This shows 
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the key role that coupling between q+G and q+G' wave 
vectors with G 7^ G' plays in the plasmon-decay mecha- 
nism in Si , which is a consequence of the presence of 
directional covalent bondings in this material. Consider- 
able broadening of plasmon peaks in Si by the effect of 
local-field effects has been reported recently by Lee and 
Chang H]. 

Next, we calculate the energy-loss function in Al and 
Si for large wave vectors of |q -I- G| = 1.65qF and 
|q -I- G| = l.Qbqp, respectively, along the directions 
(1,5/6,5/6) and (1,1,1). The results we have obtained 
within RPA and TDLDA are plotted in Figs. 8 and 9, 
showing in the case of Al the experimentally determined 
double-peak structure which was recently explained 
1^ as originated from band-structure effects. 

We note that splitting of the band structure of Al over 
the Fermi level opens new channels for the creation of 
electron-hole pairs, which results in an energy-loss func- 
tion that is, at low frequencies, higher than in the case of 
a FEG (see Fig. 8). In the case of Si, the presence of the 
band gap reduces the energy-loss function (see Fig. 9). 
As a consequence, both the absorption probability and 
the electronic stopping power of Si will be, at low en- 
ergies, smaller than those of Al, though they both have 
about the same valence electron density. 



IV. ELECTRONIC STOPPING POWER 

The electronic stopping power of a moving particle of 
charge Zie is simply the energy that it looses per unit 
path-length due to electronic excitations in the solid. 
Hence, one writes 



(32) 



where Pq represents the probability per unit time for the 
projectile to transfer momentum ?iq to the electron gas, 
the sum is extended over all momentum transfers, v is 
the projectile velocity, and Hlu is the energy transfer. As 
long as recoil can be neglected, huj — Tiq- v. 

In the case of a homogeneous FEG the probability Pq 
is given by En. (p4[), and introduction of this expres- 
sion into Eq. (^2|) yields the well-known formula for the 
stopping power of a FEG Q|. For periodic crystals, the 
energy loss depends on whether the projectile is moving 
at a definite impact parameter or not. For random tra- 
jectories, introduction of Eq. (|25| ) into Eq. (|3^ ) results 
in the so-called random stopping power: 



dE 
dx 



random 



I dq^ 

BZ G 



xlm — e 



(33) 



where the integral runs over q vectors within the first 
Brillouin zone, and cj = (q + G) • v. 



In the case of charged particles moving with constant 
velocity v on a definite trajectory at a given impact pa- 
rameter b, the probability per unit time for the projec- 
tile to transfer momentum ?i(q + G) to the electron gas 
is given b y E g. (p^), and introduction of this expression 
into Eq. (|3^ ) yields the following result for the position- 
dependent stopping power: 



dE 
dx 



X-Wq+G-HKlm \-(-c^,G+K.{^-,^) 



(34) 



where the sum is restricted, as in Eq. (p9[), to those 
reciprocal-lattice vectors that are perpendicular to the 
projectile velocity, and w = (q -I- G) • v, as in Eq. (|3^). 

The main ingredient in the calculation of both ran- 
dom and position-dependent stopping powers is the in- 
verse dielectric matrix eQ^Q, (q, w), which has been dis- 
cussed in the previous section. As the maximum energy 
that the moving particle may transfer to the target is 
{TiLj)max — Qv, the number of bands that are required 
in the evaluation of the polarizability Xg G'(qi ^) '^^ 
(^) depends on the projectile velocity. Well-converged 
results have been found for all projectile velocities under 
study, with the use in Eq. (^) of 60 bands for Al and 200 
bands for Si. The sums over reciprocal-lattice vectors 
in Eqs. ( p3| ) and (^^ have been extended over 15 val- 
ues of reciprocal-lattice G vectors, the magnitude of the 
maximum momentum transfer being 2.9qF and 2.1qp for 
Al and Si, respectively. All calculations have been per- 
formed with full inclusion of crystalline local-field effects, 
i.e., by inversion of the full dielectric matrix eGGi{ci,uj), 
and contributions from these so-called local-field effects 
to the random stopping power of Al and Si have been 
found to be within 0.5% and 1%, respectively. 



A. Random stopping power 

In Fig. 10 we show, as a function of the projectile ve- 
locity, our full RPA and TDLDA results for the random 
stopping power of valence electrons in real Al for pro- 
tons {Zi = 1), together with the corresponding result for 
the stopping power of a FEG with an electron-density 
parameter equal to that of Al (r^ = 2.07) These 
results have been found to be insensitive to the choice of 
the projectile-velocity direction. 

As the energy-loss function of real Al is, at low frequen- 
cies, slightly enhanced with respect to the corresponding 
FEG calculation (see Fig. 8), the stopping power of the 
real target is, for projectile velocities smaller than the 
Fermi velocity and within both RPA and TDLDA, higher 
than that of a FEG by about 7%. 

At velocities over the plasmon-threshold velocity for 
which plasmon excitation becomes possible, contribu- 
tions to the stopping power come from both plasmon 
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and electron-hole excitations. These contributions have 
been calculated separately [Q, showing that contribu- 
tions from losses to plasmon excitation are independent 
of the detailed band structure of the crystal. As for the 
contribution from the excitation of electron- hole pairs, 
band-structure effects in Al are found to lower the stop- 
ping power of electrons in a FEG by about 10% at and 
just above the plasnion-threshold velocity. 

At high velocities, well above the stop ping maximum, 
the sum over the frequency oj in Eq. ( |33|) can be replaced 
by an integration from to c», and the sum rule of Eq. 
( ^ ) results in a stopping power which depends on the 
average electron density no and not on the details of the 
band structure of the target material: 



dE 
dx 



random 



no In ■ 



(35) 



While at low velocities the contribution to the total 
energy loss due to excitation of inner-shell electrons is 
negligible small, at velocities larger than the Fermi ve- 
locity it is necessary to allow for this contribution. The 
cross sections for the ionization of inner shells in Al were 
obtained by Ashley et al in the first Born approx- 
imation utilizing atomic generalized oscillator strength 
functions. By adding the contribution from core elec- 
trons to that of valence electrons (this contribution does 
not depend, at these velocities, on the details of the band 
structure) these authors found a nice agreement with ex- 
periment. Good agreement with experiment was also 
shown in Ref. ||] by adding the energy loss from core 
electrons in Si, as taken from Walske's calculations [ p5[ , 
to that from valence electrons. 

Our full calculation of RPA and TDLDA stopping pow- 
ers of valence electrons in Si is plotted in Fig. 11, as a 
function of the projectile velocity, with Zi ~ 1, together 
with the corresponding result for the stopping power of 
a FEG with rg — 2.01 53|. As in the case of Al, these 
results have been found to be insensitive to the choice of 
the direction of the projectile velocity. 

The band gap in Si clearly makes the stopping power 
of this material different from that of a simple metal as 
Al. Though the average band gap of Si is small compared 
to the bandwidth of both valence and conduction bands, 
there are fewer low-energy excitation levels available than 
in the case of a FEG, and a lower energy loss is expected. 
Our calculations show that at low velocities the stopping 
power of Si is smaller than that of a FEG with rs — 2.01 
by about 10%. Also, the stopping power of Si is found not 
to be, at very small projectile velocities, proportional to 
the velocity, in agreement with experimental low-velocity 
stopping powers of this material 1 5^ . 

As in the case of Al, contributions from plasmon ex- 
citation to the stopping power of Si are found to be in- 
sensitive to the band structure, while contributions from 
electron-hole excitation yield a stopping power of this 
material which is just over the plasmon threshold about 
20% lower than that of a FEG. At high velocities, well 



above the stopping maximum, all calculations are found 
to converge. 



B. Position-dependent stopping power 

We have carried out, from Eq. (^), calculations of 
the position-dependent electronic energy loss of protons 
in Al and Si. In the case of Al, calculations have been per- 
formed for slow ions moving along the (100) and (111) di- 
rections , showing that the existence of small electron- 
density variations in this material result, through the 
presence of non-negligible off-diagonal elements in the in- 
teracting density-response matrix, in position-dependent 
corrections to the random stopping power of up to 10% 
and 20%, respectively. The maxima in the stopping 
power for trajectories along the interstitial regions and 
the minima near the cores are associated with corre- 
sponding maxima and minima in the integrated elec- 
tronic densities along the projectile trajectories. Never- 
theless, a local-density approximation (LDA), according 
to which the position-dependent stopping power is ob- 
tained as that of a homogeneous electron gas with an 
electron density equal to the average electron density 
along the projectile trajectory, is found to predict cor- 
rections to the random stopping power which are, for 
low projectile velocities, too small. 

A contour density plot of the square lattice containing 
both the projectile trajectory along the (110) direction 
in Si and the impact vector b is displayed in Fig. 12. 
The integrated density in the (110) channel varies from 
rg = 1.49 at the atomic row (full diagonal of the fig- 
ure) to rg = 3.37 at the center of the channel (solid line 
with an arrow) , thereby showing valence-electron density 
variations of up to 80% typical of covalent crystals like 
silicon. 

In Fig. 13 we show, as a function of the projectile 
velocity, our full RPA calculation of the stopping power 
of Si for best channeled ions moving along the (110) di- 
rection. The stopping power of valence electrons in Si 
for these channeled ions is found to be, at intermediate 
velocities, about 20% smaller than the random stopping 
power. On the other hand, we note that the stopping 
maximum for channeled ions is located at the same ve- 
locity as in the case of random ions, while in the case of a 
homogeneous electron gas with an electron density equal 
to the average electron density along the channel would 
be located at a lower value of the velocity, as shown by 
the LDA calculation represented by a dotted line. We 
also note that within the LDA position-dependent cor- 
rections to the random stopping power for slow ions are 
predicted to be too small; for velocities above the stop- 
ping maximum, the local-density approximation yields 
unrealistic values for the energy loss of best channeled 
ions, which are far below our full band-structure calcula- 
tions. At high velocities, position-dependent and random 
stopping powers are found to converge. 
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V. CONCLUSIONS 

We have presented a survey of current investigations of 
the impact of band structure effects on plasmon energies, 
plasmon hnewidths, dynamic structure factors, and both 
random and position-dependent stopping powers of Al 
and Si. 

New ab initio calculations for the plasmon lifetime in 
Al and Si have been presented. In the case of Al, we 
find the full width at half maximum to be in the limit as 
q ^ of 0.2cV, below the early calculation carried out 
by Sturm with use of the nearly free electron approxima- 
tion 1^. The plasmon linewidth is found to increase as 

for small values of the momentum transfer q, while it 
increases very quickly near the plasmon cutoff, in agree- 
ment with the experiment. A similar behaviour is found 
in the case of Si, though plasmon linewidths in this ma- 
terial are found to be about five times wider than in Al. 

As for the random stopping power, we find that it is, 
at low velocities, smaller in Si than in Al, though they 
both have nearly the same valence electron density. At 
velocities just over the plasmon threshold they are both 
below the stopping power of a FEG with ~ 2. At 
high velocities there are no band structure effects. The 
random stopping power of Si is found to be, at low and 
intermediate velocities, about 10% smaller than that of 
Al, in agreement with experimental measurements for ei- 
ther protons 1^ or antiprotons |Q . 

Differences between random and position-dependent 
stopping powers of Al are found to be up to 10% for pro- 
jectiles incident in the (100) direction and up to 20% for 
projectiles moving in the (111) direction. As for Si, the 
stopping power for best channeled ions along the (110) 
direction is found to be diminished with respect to the 
random stopping power by about 20% at velocities near 
the stopping maximum. 
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FIG. 1. The RPA energy-loss function for Al, for several 
wave vectors along the (100) direction: q — (0.2, 0, 0)27r/a 
(solid line), q = (0.4, 0, 0)27r/a 

(long-dashed line), q = (0.6, 0, 0)27r/a (short-dashed line), 
and q — (0.8, 0, 0)27r/a (dashed-dotted line). The vertical 
dotted lines represent the corresponding results for a FEG 
with ra = 2.07, which are simply S functions. 

FIG. 2. RPA (solid line) and TDLDA (dotted hue) plas- 
mon excitation energies for wave vectors along the (100) di- 
rection in Al. The corresponding calculations for a FEG are 
represented by dashed and dashed-dotted lines, as obtained 
within RPA and TDLDA, respectively. The experimental re- 
sults are represented by triangles 

FIG. 3. RPA plasmon linewidths of Al for wave vectors 
along the (100) direction. Solid circles and open triangles 
represent calculated linewidths without and with inclusion of 
plasmon decay through excitation of two electron-hole pairs. 
Solid triangles represent the experimental results. 

FIG. 4. Scaled RPA (open triangles) and TDLDA (open 
circles) plasmon linewidths of Al for wave vectors along the 
(100) direction, as obtained with inclusion of plasmon decay 
through excitation of two electron-hole pairs. Solid triangles 
represent the experimental results. 

FIG. 5. The RPA energy-loss function for Si, for several 
wave vectors along the (111) direction: q = (0.2, 0.2, 0.2)27r/a 
(solid line), q = (0.4, 0.4, 0.4)27r/a (dashed line), and 
q — (0.6, 0.6, 0.6)27r/a (dashed-dotted line). The vertical lines 
represent the corresponding results for a FEG with rs = 2.01, 
which are simply 5 functions. 

FIG. 6. RPA (solid line) and TDLDA (dotted hue) plas- 
mon excitation energies for wave vectors along the (111) di- 
rection in Si. The corresponding RPA calculations for a FEG 
are represented by dashed lines. The experimental results are 
represented by triangles. 

FIG. 7. RPA (solid circles) and TDLDA (open triangles) 
plasmon linewidths of Si for wave vectors along the (111) di- 
rection. Crystalline local-field effects have not been included 
in this calculation. 



FIG. 8. The RPA (solid lines) and TDLDA (dashed hues) 
energy-loss function for Al, for a large wave vector of magni- 
tude |q-|- G| = 1.65gF along the direction (1,5/6,5/6). Plain 
solid and dashed lines represent RPA and TDLDA energy-loss 
functions of a FEG with = 2.07. 
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FIG. 9. The RPA (solid lines) and TDLDA (dashed lines) 
energy-loss function for Si, for a large wave vector of mag- 
nitude |q + G| = 1.95qF along the direction (1,1,1). Plain 
solid and dashed lines represent RPA and TDLDA energy-loss 
functions of a PEG with rs = 2.01. 



FIG. 10. Full band-structure calculation of the random 
stopping power of valence electrons in Al for protons {Zi — 1), 
as a function of the projectile velocity. Solid and open circles 
represent the results obtained in the RPA and the TDLDA, 
respectively. RPA and TDLDA stopping powers of electrons 
in a FEG with rs = 2.07 are represented by solid and dashed 
lines, respectively. 

FIG. 11. Full band-structure calculation of the random 

stopping power of valence electrons in Si for protons [Zi — 1), 
as a function of the projectile velocity. Solid and open circles 
represent the results obtained in the RPA and the TDLDA, 
respectively. RPA and TDLDA stopping powers of electrons 
in a FEG with rs = 2.01 are represented by solid and dashed 
lines, respectively. 

FIG. 12. Contour density plot of the valence electron den- 
sity in the plane defined by the (100) and (010) vectors in 
Si. 

FIG. 13. Full band-structure calculation of the RPA stop- 
ping power of valence electrons in Si for best channeled pro- 
tons (Zi = 1) in the (110) (solid circles), as a function of the 
projectile velocity. Solid and dashed lines represent the ran- 
dom stopping power of real Si and of a FEG with Vs = 2.01, 
respectively. The dotted line represents the stopping power 
of a FEG with = 3.37. 



10 



50 



40 



30 



Figure 1 



T I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



3 



20 



10 







I- j- 4 ^ .-I 




10 




-I -^1 - I ~ I - I- I I U 



-I I -- j-^ 



15 20 

co(eV) 



25 



Figure 3 



T I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



> 



4 



o 



in 



-V 



q 







J I I I I I I I I I I I I I I \ I I I I I I I I I L 



c 

J L 



0,1 0,2 0,3 0,4 0,5 0,6 



Figure 4 

T I I I I I I I I I I I I I I I I I I I I I I I I I I r 




I I I I I I I I I I I I I I I \ I I I I I I I I I \ I I I 

0,1 0,2 0,3 0,4 0,5 0,6 



Figure 5 




10 15 20 25 30 

co(eV) 



Figure 6 




0,2 0,4 0,6 0,8 1 

q(27r/a) (1,1,1) 



3,5 



> 

I ^ 

o1,5 

E 1 



Figure 7 



0,5 



T I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



J I I I I I I I I I I I I I I \ I I I I I I I I I L 



J L 



0,1 0,2 0,3 0,4 0,5 0,6 



Figure 8 




20 40 60 80 100 

co(eV) 



Figure 9 

0,15 I I I I I I I I I I I I I I I I 




20 40 60 80 100 120 

co(eV) 



Figure 10 




0.5 1 1.5 2 2.5 3 3.5 4 

V (a.u.) 



Figure 11 




0,5 1 1,5 2 2,5 3 3,5 4 

V (a.u.) 



Figure 12 




Figure 13 




0,5 1 1,5 2 2,5 3 3,5 4 

V (a.u.) 



